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Abstract 

Spin precession was nonthermally induced by an ultrashort laser pulse in orthoferrite DyFeC>3 
with a pump-probe technique. Both circularly and linearly polarized pulses led to spin precessions; 
these phenomena are interpreted as the inverse Faraday effect and the inverse Cotton-Mouton ef- 
fect, respectively. For both cases, the same mode of spin precession was excited; the precession 
frequencies and polarization were the same, but the phases of oscillations were different. We have 
shown theoretically and experimentally that the analysis of phases can distinguish between these 
two mechanisms. We have demonstrated experimentally that in the visible region, the inverse Fara- 
day effect was dominant, whereas the inverse Cotton-Mouton effect became relatively prominent 
in the near-infrared region. 

PACS numbers: 78.20.Ls, 75.50.Ee, 75.40.Gb, 78.47.J-, 
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I. INTRODUCTION 



Magnetization switching triggered by femtosecond laser pulses has been studied in re- 
cent years. Ultrafast demagnetization in ferromagnetic metals and semiconductors has also 
been reported.-^ These phenomena show thermal magnetic switching with light pulses on 
picosecond time scales.- However, heat-assisted spin reorientation is relatively slow because 
of the thermal diffusion time. 

A light pulse with a certain polarization nonthermally modifies the electron spin state.-^ 
Recently, it has been reported that spin precession is induced by a circularly polarized pulse 
in antiferro magnetic (AFM) DyFeC>3 with weak ferromagnetic (FM) moment.- The phase 
of spin precession changes by 180° on reversal of the pump helicity. The interpretation of 
this phenomenon is that an effective magnetic field pulse parallel to the pump wave vector is 
induced by the circularly polarized light pulse, giving rise to the precession. The magnetic 
field generation effect is referred to as the inverse Faraday effect (IFE). The same effect has 
also been observed even in pure AFM NiO with no net magnetic moment in the ground 
state.- The resonance frequencies of AFM materials reach the terahertz range, which is 
several orders of magnitude higher than that of FM materials. For that reason, AFM 
materials attract much attention in the context of ultrafast spin control.- - — Spin precession 
is also observed with a linearly polarized pump pulse, in particular, a pulse polarized in a 
direction non-parallel to the crystal axes. This phenomenon is called the inverse Cotton- 
Mouton effect (ICME).— A detailed review of these phenomena can be found in Ref. 

The ultrafast IFE and ICME are interpreted as impulsive stimulated Raman scattering 
(ISRS).— ~— An electron in the ground state is excited by the pump pulse into a virtual state, 
which changes the orbital momentum of the electron. The nonzero orbital momentum flips 
the electron spin with spin-orbit coupling in the virtual state. The excited electron radiates 
a photon and transits to the final state. The energy gap between the final and ground states 
corresponds to the spin precession energy. 

ISRS is a modulation of the dielectric permittivity by the pump pulse and should be 
dependent on the properties of the pulse, such as its polarization, wavelength, and fluence. 
Therefore, examining the dependence of the photoinduced spin precession on these proper- 
ties will help us to understand the ISRS mechanism. In particular, it is not obvious how 
the pump photon energy influences spin precession. An action spectrum of photoinduced 
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spin precession should indicate the relation between the optical excited state and the spin 
precession via ISRS. 

In the majority of previous publications, the excitation of spin oscillations by ultrashort 
laser pulses was associated with IFE and ICME separately. In the present work, we report 
spin precession induced via ISRS as functions of the pump pulse polarization and wavelength. 
We found that both effects, IFE and ICME are working in the same way, exciting the same 
mode of spin precession. The phases of the spin precession via IFE and ICME differ by 
90°, allowing the two effects to be distinguished. We found an essential dependence of the 
phase on the pump wavelength and demonstrated that the IFE and ICME are dominating 
effects in different spectral regions, in the visible region and in the near-infrared region, 
respectively. Thus, the analysis of the phase difference of the spin precession reveals the 
mechanism of ISRS. 



II. PHYSICAL PROPERTIES 

A. Crystallographic and magnetic properties 

DyFe03 is a rare-earth orthoferrite and crystallizes in an orthorhombic structure 
(Pbnm).— Spins of the Dy 3+ ions are not ordered above 4 K. Four Fe 3+ ions occupy positions 
(1/2, 0, 0), (1/2, 0, 1/2), (0, 1/2, 1/2), and (0, 1/2, 0) in the unit cell. In the exchange ap- 
proximation, the arrangement of their magnetic moments, M\, M 2 , M 3 , -M 4 , corresponds 
to one of the four patterns G { : M u = -M 2i = M 3i = -M iu F 4 : M u = M 2i = M 3i = M 4i , 
Ai : M u = -M 2i = -M 3i = M 4i , and Q : M u = M 2i = -M 3i = -M Ai {% = x, y, z). 
DyFeOa crystal has the spin arrangement r4(G x A y F z ) and belongs to the magnetic point 
group m'm'm above the Morin point and below the Neel temperature, at 37 K < T < 
Tjy = 645 K.— ~— Because of the superexchange interaction, the spins are almost completely 
arranged antiferromagnetically along the x-axis. Due to the Dzyaloshinskii-Moriya interac- 
tion, all spins tilt by about 0.5° toward the z-axis.— >22 Usually the conditions Mi ~ M 3 
and M 2 ~ M4 are valid and a simpler model with just two different sublattice magnetic 
moments, Mi and M 2 , with \Mx\ = \M 2 \ = M , can be employed.—^ In what follows, 
this two-sublattice model will be used. We denote the FM vector by M = Mi + M 2 and 
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the AFM vector by L = Mi — M^- These vectors are subject to constraints 



(M ■ L) = 0, M 2 + L 2 = 4M 2 . 



(1) 



The dynamics of M(t) and L{t) is described by Landau-Lifshitz equations 



26.31.32 



dM(t) 

dt 

dL(t) 
dt 



7 {[M(t)x/ i ctt ] + [ J L(t)xi? ctt ]}, 



7{[M(t) x H ca ] + [L(t) x h etl }}, 



(2) 



(3) 



where 7 = g^/fr (> 0) is the gyromagnetic constant, /ib is the modulus of the Bohr 
magneton, g is the gyromagnetic ratio, g « 2 for orthoferrites, and i2" cff and fo efr are the 
effective magnetic fields. Using the magnetic energy of an orthoferrite, the effective fields are 
denoted as H cR = —dT-L/dM and h cS = —dT-L/dL, where the Hamiltonian is given by2&^ 



The last term describes the Dzyaloshinskii-Moriya interaction and d is parallel to the y-axis. 
Equations (j2J), (j3J), linearized above the ground state determined by the energy (TJ]), yield 
eigenmodes of oscillations of the vectors M(t) and L(t). These spin precession modes for 
the T 4 ground state with the equilibrium values of M z 7^ and L x 7^ (see Fig. Ha)) are 
described as follows 
with as following: 



where z and x are unit vectors parallel to the z-axis and x-axis, respectively, and the 
variables m(t) and l(t) correspond to two eigenfrequency modes, as shown in Fig. [TJ[b) and 
(c). The components m x , m y and l z oscillate at the quasi-ferromagnetic resonance (F-mode) 
with the angular frequency up. On the other hand, l x , l y , and m z oscillate at the quasi- 
antiferro magnetic resonance (AF-mode) with the angular frequency waf-~ Those resonance 
frequencies are given by^— 



(4) 



M(t) = M z z + m(t), 



L{t) = L x Z+l(t). 



(5) 
(6) 



co F = 2'yM oy /A(q 3 -q 1 ), 



(7) 
(8) 



u AF = 2 7 M 0V / rf 2 -^i, 
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where Mq = \M\\ = I-M2I, and the anisotropy constants p\ and pz are omitted, because 
their contribution to frequencies is negligible. (It is worth to note here, that the exchange- 
relativistic constant of the Dzyaloshinskii-Moriya interaction has the same order of magni- 
tude as a square root of the product of exchange and relativistic constants like and 
thus should be kept in the above expressions.) The orbits of the spin precession and the 
temporal response of M and L are different in the two modes. 

B. Optical and magneto-optical properties 

DyFeC>3 has the d-d transitions & Ai g — > A E g , 4 Ai g centered at the wavelength of 500 nm, 
6 Ai g — > A T2 9 at 700 nm, and 6 A\ g — > A T\ g at 1000 nm.S This crystal is optically biaxial, so 
the components of the dielectric permittivity tensor are e xx 7^ e° 7^ e° zz . The birefringence p 
stems from the difference in the refractive indices, such as An xy = n x —n y (p = 2irAn/X). On 
the other hand, the magnetization M z leads to the Faraday rotation for the light propagating 
along the z-axis. The Faraday effect is much smaller than the effect of birefringence. In 
DyFeOs, the Faraday rotation <fi and the birefringence per unit length are = 1.6 x 10 3 
deg/cm and p = 1.2 x 10 5 deg/cm, respectively, at 800 nm.-- 

C. Interaction of the light pulse and the medium 

The interaction of the magnetic medium and transmitting light is described by the dielec- 
tric permittivity tensor 8^.—^- By virtue of the Onsager principle, if absorption is negligible, 
£jj can be divided into antisymmetric and symmetric parts, (e^ = — and (e|- = e^), 
with real and imaginary components, respectively. For a transparent medium, the tensor 
components can be written in the following general form 

£ ij = £ ij + ifijkMk + ig ijk L k + a m M k Mi + b ijkl L k Lt + c ijM M k L h (9) 

where is a magnetization-independent term having a symmetric part only. By taking into 
account the symmetry of orthoferrite, the terms in the first line (except the ejj ) represent 
the antisymmetric part of Eij, and the terms in the second line describe the spin-dependent 
symmetric part of the permittivity tensor. The symmetry of the fourth rank tensors a^z, 
bij k i, and c^u is determined by the magnetic and crystal point groups, and fij k and g^ k are 
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the third rank tensors, antisymmetric over the first pair of indices, e.g., fijk = —fjik- Tensors 
CijM and gijk originate from the Dzyaloshinskii-Moriya interaction. The Hamiltonian of the 
interaction between the light pulse and the medium in SI unit isM- 

n mt = £ fm£*(t), (io) 

where £j(t) is the time-dependent amplitude of the light in the pulse. A circularly polarized 
pulse propagating along the z-axis can be described in the form (£ x (t),£ y (t)) = t^(1, 
where the ± indicate the opposite senses of helicity. A linearly polarized pulse with the 
polarization inclined on at an angle 9 with respect to the x-axis can be described in the form 
(S x (t), £ y (t)) = £(t)(cos6,sin8). Then a straightforward calculation gives the Hamiltonian 
of the interaction with the medium of the form 

U£ = l£(t)£*(t) (e s xx + e s yy T 2ie a xy ) , (11) 

U V Z = \£(t)£*(t) (e s xx cos 2 6 + e s yy sm 2 6 + e s xy sm26) , (12) 

for circularly and linearly polarized pulses, respectively. Nonzero components of the tensors 
e% and e% are listed in Table UJiiii 

When a pulse is incident on a medium, the interaction between the pulse and the medium 
is given by Eqs. (IT0" ]) - (fI2]) . The incident pump pulse generates effective pulsed fields H cS = 
—dUint/dM and h cS = —dH mt /dL. Both effective fields are proportional to the intensity 
of the light, £{£)£* (t). If the pulse duration A is much shorter than the period of spin 
oscillations, A I/cjf, 1/^af, the real pulse shape can be replaced by the Dirac delta 
function, £{t)£*{t) — > Io5(t), where Io = J £{t)£*{t)dt is the integrated pulse intensity. The 
light-induced effective fields H cS and h cS can be regarded as being proportional to the delta 
function 5(t) as well. For a light pulse propagating along the z-axis, H cS and h eS generated 
by a circularly polarized pulse are 

H — \2*{flxxzz &yyzz)-M z -\- {c xxzx ~\- C yyzx )L x i *2f xyz ^Z, (^3) 

o 

h ~ [^ibxxxx b yyxx ) L x (cxxzx c yyzx }]\/I z i 1g xyx pi.^ (1^) 

o 

respectively. The phenomenon of generating these effective magnetic fields is known as IFE. 

For a magnetic field pulse of a short duration, the action of the light-induced effective fields 
within the delta-function approximation can be described as an instantaneous deviation of 

6 



TABLE I. The dielectric permittivity tensor e^- as a function of magnetic component. Modulation 
of the dielectric permittivity in F-mode and AF-mode is shown in columns 3 and 4, respectively. 
No external magnetic field is present, ef- = —e% and = e% 



Tensor 


Static 


F-mode 


AF-mode 


element 


(m(t)=0, i(t)=0) 


(m x ,m y ,l z ^ 0) 


(l x Jy,m z / 0) 


£ S 


&XXZZ-M z -\- b xxxx L x 
C XXZX M Z L X 





y^&xxzzMz -\- C xxzx L x )fn z 
i.'^b xxxx L x ~\~ C xxzx M z ^jl x 


£ s 

yy 


&yyzzM z byyxxLg 
~i~ C-yyzx-M Z L X 





{^CLyyZzMz; ~\- Cyyz X L x ^TTlz 
+ (2byy XX L x + CyyzX-Mz)l X 


^zz 


dzzzzM^ + bzzxxL^ 
Czzzx-M Z L X 





ficizzzzMz + c zzzx L x )m z 
{2b zzxx L x + c zzzx M z )l x 


£ S 

xy 








(2b X y X yL X + C X y Z yM Z ^ly 







(^Q'ZxxzM z + C zxxx L x )m x 
~\~ i^b zxxz L x ~\~ C ZXZZ A4 ~z)l z 





£ S 

yz 





(2(ly z y z M Z + CyzyX-Lx^Hly 





e a 
xy 


ifxyzM z + i(jxyxL x 





ifxyzlTbz ~\~ iQxyxlx 


e a 

^zx 





ifzxy m y 





e a 
yz 





i fyzx^Tlx ~i~ idyzz^z 






the FM and AFM vectors, AM = M(t = +0) - M(t = -0) and AL = L{t = +0) - 
Lit = — 0), from their equilibrium positions, M[t = — 0) = M z z and L(t = — 0) = L x -x, 
respectively. After vanishing of the pulsed effective field, the spins precess around the 
effective fields corresponding to their equilibrium directions following the Landau-Lifshitz 
equations, based on the Hamiltonian (J3J). Thus the action of the pulse can be regarded as a 
creation of some (non-equilibrium) initial conditions for the Landau-Lifshitz equations. The 
deviation of the FM and AFM vectors induced by the circularly polarized pulse is described 
by 

AM CT± = 0, (15) 
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\2(b xxxx byyxx &XXZZ Q"yy Z z) z L x i^CxXZX ~\~ Cyy zx )M z 

o 

— {.Cxxzx Cyyzx)L x -p 2(f X y Z L x Qxyx-^z)\y ■ (16) 

Here, M is not affected by the effective field directly, whereas l y of L(t) takes nonzero 
deviations. The resonance mode with l y ^ is AF-mode. In Fig. [He), two spins M i and 
M 2 , as well as their sum and difference M and L, move toward positions 2 or 4, and the 
spins precess around their ground state directions. Because the spins have only an l y variable 
component when the effective magnetic field disappears, spin precession starts at position 2 
or 4 (see Fig. Q^c)). 

Similarly, effective magnetic fields induced by a linearly polarized pulse are 

jjeflF.Ui, = _b^)_^ 2axxzzMz + Cxxzx l x ) cos 2 9 + (2a yyzz M z + c y y ZX L x ) sin 2 0]z, (17) 

h cSM = _h^[{(2b xxxx L x + c XX zxM z ) cos 2 9 + (2b yyxx L x + c yyzx M z ) sin 2 9}k 

+ (2b xyxy L x + c xyzy M z ) sin 26* • y]. (18) 

These effective magnetic fields are induced via ICME. The deviations of the FM and AFM 
vectors created by the effective field are 

AM lin = ^ 2bxyxyL 2 + CxyzyMzLx) sin2 • Z, 

(19) 



A£ iin = lh[{(2b xxxx M z L x + CxxzxM 2 z - 2a X xzzM z L x - c xxzx L 2 x ) cos 2 9 

+ (2byyxxM z L x + Cyy ZX Ml ~ 2(Xyy zz M z Lx ~ Cyy zx L x ) SVC? 9}f 



- (2b xyxy M z L x + c xyzy M 2 z ) sin 29 ■ x]. (20) 

Here, the components m z of M and l x , l y of L are affected by the effective field. This 
precession mode is also an AF-mode, but the initial direction of the spin deviation differs 
from that for the circularly polarized pulse case. 

A pulse propagating along the x- or y- axis should trigger the spin precession with both F- 
and AF- modes. The amplitude and the phase of the precession depends on the polarization 
of the pulse. 
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III. EXPERIMENTAL RESULTS 



A. Experimental setup 

We studied photoinduced spin precession in DyFeC>3 using a pump-probe magneto-optical 
technique, as shown in Fig. [2j DyFeOs single crystals were grown by the floating-zone 
method, and the orientation of the faces were determined by back-reflection x-ray Laue 
photographs.— Faces with a width of few millimeters were mechanically polished. The 
sample thickness was 140 fim, except for a thickness dependent measurement. The sample 
was placed in a cryostat at 77 K with no external magnetic field. Optical pulses with a 
central wavelength of 790 nm, a duration of 150 fs, and a repetition rate of 1 kHz were 
emitted from an amplified Ti: sapphire laser. The beam was separated into two beams 
by a beamsplitter. One was employed as the probe beam, and the other was injected 
into an optical parametric amplifier (OPA), which converted the incident beam to signal 
and idler beams, in the wavelength ranges 1140-1580 nm and 1580-2570 nm, respectively. 
Furthermore, the signal and idler beams were frequency- doubled with a /3-BaB204 (BBO) 
crystal if necessary. Then unwanted beams were cut by color filters. The ranges of the pump 
wavelength were 600-750 nm (second harmonic of the signal pulse), 850-1100 nm (second 
harmonic of the idler pulse), and 1140-1500 nm (the signal pulse). 

Figure [3] illustrates the circular and linear polarizations employed for the pump and probe 
pulses. Circularly polarized pulses are denoted as a ± . Linearly polarized pulses, denoted LI, 
L2, L3, L4, L5, and L6, were tilted at — vr/4, 7r/4, 0, 7r/2, and from the x-axis, respectively, 
where tana = 2. 

The fluence of the pump pulse was varied from 15 to 130 mJ/cm 2 , depending on the 
wavelength. The pump pulses were focused on the sample to spot sizes of 50-100 /mi. The 
probe pulses were linearly polarized and had a pulse fluence of 1 mJ/cm 2 . The probe beam 
was vertically incident on the surface of the sample, whereas the pump beam was incident 
at the angle of 7°. The transmitted probe pulse was divided into two orthogonally polarized 
pulses by a Wollaston prism, and each pulse was detected with a Si photodiode. The ratio 
of the signals from the detectors allowed us to determine the angle of the probe polarization. 
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B. Dependence of the polarization rotation on the pump pulse polarization 

Figure H] illustrates the polarization of the propagating pulses in the medium with bire- 
fringence. For the sake of simplicity, we will discuss the picture of the light propagation 
without taking the Faraday effect into account. Pulses with circular polarization or linear 
polarization nonparallel to the crystal axis are transformed, whereas pulses with linear po- 
larization parallel to the crystal axis, corresponding to the normal modes of light in the 
media, are not. Thus, for the pulses with general linear polarization or pulses with circular 
polarization, the real and imaginary parts of E x E y , which are responsible for the terms in 
Eqs. (fI5]) - (pD]) including sin 20, and f xyz and g xyx , respectively, will oscillate in space along 
the pulse propagation direction, while they remain uniform only for pulses linearly polarized 
parallel to the crystalline axis. Therefore, the effective magnetic field and spin precession 
generated by IFE and ICME will be different at different positions in the sample. 

Figure |5] shows the polarization rotation of the probe pulse as a function of the delay 
time between the pump and probe pulses. The pump wavelength was 1050 nm, and the 
polarizations were a , LI, L2, L3, and L4. The probe polarization was L4. When the pump 
polarizations were a , LI, and L2, oscillation of the polarization rotation was observed. The 
frequency of the oscillation was 210 GHz at the temperature T = 77 K, in agreement with 
previous infrared and Raman experiments.— 1 ^ 1 ^ In Figs. 0(c) and (d), the pump pulses with 
polarizations L3 and L4 did not induce oscillation of the probe polarization. Polarizations 
a ± , LI, and L2 had E x E y components, but L3 and L4 did not, as shown in Fig. HI 



C. The influence of magnetization on the probe polarization 

Modulation of the dielectric permittivity leads to oscillation of the probe polarization in 
the sample. The origin of the probe polarization change can be attributed to the Cotton- 
Mouton effect and the Faraday effect. The Cotton-Mouton effect is magnetic linear bire- 
fringence based on , whereas the Faraday effect is magnetic circular birefringence based 
on e a xy . 

In order to identify the effect giving rise to the polarization rotation as observed in Fig. 
we set <j ± for the pump polarization and L5 and L6 for the probe polarization. For L5 and 
L6, the Faraday effect leads to rotation of the probe polarization in the same direction 
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for both probe polarizations, whereas the Cotton-Mouton effect leads to rotation in the 
opposite direction. Therefore, the dominance of the rotation of the probe polarization can 
be distinguished. Figure shows that the polarization rotations of two probe pulses with 
polarizations L5 and L6 oscillated in the same direction. This indicates that the contribution 
of the Faraday effect is dominant and that of the Cotton-Mouton effect is negligible for the 
probe wavelength of 800 nm. This is consistent with the fact that IFE is dominant for the 
pump wavelength of 800 nm (see below). 



D. Dependence of the polarization rotation on the pump wavelength 

The oscillation of the probe polarization originates from spin precession. Therefore, 
the phase of the oscillation indicates the direction of an effective magnetic field induced 
by the pump beam. The dependence of the effective magnetic field and reorientation of 
magnetization on the pump wavelength gives information about the interaction of the light 
pulse and the magnetic medium. 

An experiment was performed with four types of pump polarizations, a ± , LI, and L2. 
The differences between the oscillations for a + and a~ and between those for LI and L2 were 
measured. Figure [7] (a) shows the initial phase £ of the oscillation of the probe polarization 
versus pump wavelength. The oscillation is described by 9(t) = Asin(ujt + £) at t > 0, where 
A is the amplitude, u is the angular frequency, and £ is the initial phase. The initial phase 
was close to 0° (or ±180°), when the pump wavelength was 800 nm. This is consistent with 



Ref. 



45 



When the pump wavelength was between 1000 nm and 1100 nm, the initial phase was 
closer to ±90°. When the pump wavelength was above 1200 nm, the initial phase was be- 
tween 0° and 90°. By comparing two samples with thicknesses of 140 /im and 170 jum, it 
was confirmed that the sample thickness does not affect the phase shift (data not shown). 
Figure [7J (b) represents the amplitude A of the oscillation as a function of the pump wave- 
length. The amplitude A is proportional to the pump fluence, thus justifying normalization 
of the amplitude by the fluence. Because of the transformation of the pulse polarization in 
the medium, the magnetic field differs at different positions. Therefore, the amplitude A 
was not simply proportional to the magnitude of the generated magnetic field. However, 
when the pump wavelength was from 700 nm ( 6 Ai — > 4 T 2 ) to 1000 nm ( 6 Ai — > 4 Tx), the 
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amplitude was larger than that of the other region in Fig. [7] (a) . This result suggests that 
the photoinduced spin precession is related to the electron transition. 



E. Pump probe measurement in (100) and (010) oriented crystals 

To determine all dielectric permittivities, we performed pump-probe measurements in 
(100) and (010) oriented crystals. The pump wavelength was 750 nm, and the crystal 
thickness in both cases was 100 /im. However, in contrast to the previous experiments,- 
oscillation of the polarization of neither F- nor AF-modes was observed in either propagation 
direction. 



F. The dependence of polarization rotation on temperature 

It is well known that magnon frequencies in orthoferrites strongly depend on the 
temperature .-^^^ We measured the temperature dependence of the spin precession prop- 
erties in DyFeC>3. The frequencies of the oscillations for pump wavelengths of 750 nm and 
1200 nm are shown in Fig. [8] (a), in comparison with previously reported spin precession.— 



Our data show excellent agreement with Refs. If] and |31|, regardless of the pump wavelength. 
The frequency decreases with approaching the Morin point T r = 37 K because of magnon 
softening associated with the spin reorientation.— The temperature dependence of the 
initial phase £ of the spin precession for pump wavelength of 750 nm is shown in Fig. [8] 
(b). The initial phase was close to 0° or 180° with a jump at T=150 K. It is worth to note 
that at this temperature the frequencies of F-mode and AF-mode become equal, that is 
A(q 3 — qi) = d 2 — Aq\. Furthermore, the energies of two domain walls with the spin rotation 
in (010) and (001) planes become equal at this point, which leads to the reconstruction 
of domain walls.— However, we were not able to find the relation between the properties 
described above and the initial phase shift. 
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IV. DISCUSSION 



A. Landau Lifshitz equations 

According to the results of the previous section, the number of essential dielectric per- 
mittivity components can be reduced. First, we found that pump-probe measurement in 
(100) and (010) oriented crystals revealed that e s zz , e s xz , e yz , e% z , and e yz were negligible. In 
addition, pump pulses with L3 and L4 polarizations did not trigger spin precession in the 
(001) oriented crystal in Fig. |5j Polarizations L3 and L4 had only electric field components 
E x and E y , respectively. Therefore, the terms containing cos 2 9 and sin 2 6 in Eqs. (IT7|) . f fTBj) 
and ( 12 Op were also negligible. 

Moreover, it has been reported that f xyz M z and g X yxL x are of the same order of magnitude 
for orthoferrites.— In contrast to that, for the AF-mode the ratio of m z and l x is \m z /l x \ = 
\L X /M Z \ ~ 100. Thus, f xyz m z 3> g xyx l x , and one can ignore the term g xyx l x . In addition, 
Fig. |6] indicates that the observed oscillation of the polarization was dominated by the 
imaginary part of the dielectric permittivity e xy = if xyz (M z + m z ) + ig xyx {L x + l x ). Because 
fxyz m z 3> g xyx l x , the phase of m z corresponds mostly to that of e a xy and that of the oscillation 
of the polarization. 

These findings simplify the dielectric permittivity tensor. In Table [U the tensor elements 
in the AF-mode column are proportional to l x and negligible, except for e s = (b xyxy L x + 
c xyzy M z )l y and e xy = if xyz m z . Here we suppose that a pulse is incident on a (001) oriented 
crystal. We can simplify the effective magnetic field and the dynamics of the magnetization 
induced by the circular polarization: 

H eS ^ = T W 4 }/ ^ z, (21) 

h cS ' a± = 0, (22) 

AM a± = 0, (23) 

AL a± = ± 2Mp^: y . (24) 
In the case of the linear polarization one in turn obtains: 

H cSM = 0, (25) 

& eff,lin = - I -^-(2b xyxy L x + c xyzy M z ) sin2# • y, (26) 
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AM iin = J°( 2bxyxyL 2 + CxyzyMz L x ) sin 29 • z, (27) 

Ai lin = -lh{2b xyxy M z L x + c xyzy M 2 z ) sin 29 ■ x. (28) 

The second terms are much smaller than the first ones in Eqs. (126|) . (1271) and (1281) . re- 
spectively. As a result, IFE and ICME are induced by the contributions of e xy and e s xy , 
respectively. 

Equation ff24^) indicates that the circular polarization causes the AFM component l y and 
rotation torque of the AF-mode. On the other hand, the FM component does not change 
in Eq. ( |23l) . As a result, IFE leads to oscillations proportional to sina^F^- On the other 
hand, Eqs. ( 1271) and (1281) indicate that linear polarization causes components m z and l x . As 
shown in Fig. [TJ m z and Z x have the same phase, so ICME leads to oscillations proportional 
to coswaf^- Therefore, the initial phases of m z excited by IFE and ICME differ by 90°. 

Because the phase of m z is nearly equal to that of the oscillation of the polarization, we 
can estimate the phase of spin precession from the result in Fig. [7J Since the polarization 
of the pump pulse is transformed by birefringence, the effective magnetic field and spin 
precession differ at different positions in the medium, as shown in Fig. HI However, if one of 
IFE and ICME is dominant and the other is negligible, the time dependence of m z and the 
oscillation of the probe polarization are proportional to sin cuaf^ or cos cuaf^, respectively. 



B. Sigma model 

Nonlinear sigma model is a convenient tool for the description of linear and especially 



26 for details. It is based on the 



non-linear spin dynamics of antiferromagnets, see Ref. 
dynamical equation for the vector L only that is of the second order in time derivatives, 
whereas the vector M is a slave variable so it can be expressed through the vector L and 
its time derivative. Recently, two alternative scenarios of laser-induced excitations of spin 
oscillations in antiferromagnets have been discussed within the framework of this model. 
The so-called inertial mechanism has been proposed for canted antiferromagnets and has 
been realized experimentally for holmium orthoferrite.— Within the sigma-model approach, 
the inertial mechanism is associated with an action of the laser-induced pulse of the magnetic 
field on the vector L as a pulse of force on the massive particle. In this mechanism, the laser 
pulse creates an initial value of the time derivative, dL/dt, that in principle can lead to quite 

11 



large deviations of the vector L after the action of the pulse. In the alternative mechanism, 
the time derivative of the effective magnetic field plays a role of the driving force, leading to 
an initial deviation of the vector L from its equilibrium direction.-^ For this field- derivative 
mechanism, the amplitudes of spin deviations are expected to be smaller than for inertial 
mechanism, but can be realized for any antiferro magnet, even a purely compensated one. 
The latter mechanism has beed observed experimentally in AFM nickel oxide, where the 
Dzyaloshinskii-Moriya interaction is forbidden by symmetry- It is interesting to understand 
which mechanism is responsible for the spin oscillations observed in the present work. 
The Lagrangian density of the sigma model can be written as follows:— 

c - WamI (§ ) 2 + ^ImI { H ■ (f x £ )) + > ' < £ x d » " w "< £ ) ■ < 29 > 

where H is the effective magnetic field and W a (L) is the effective anisotropy energy 
that includes L-dependent terms from the Hamiltonian (J4]) and a contribution from the 
Dzyaloshinskii-Moriya interaction, see Eq. fl3T|) below. The slave variable, magnetic mo- 
ment M , can be easily expressed via vector L and its time derivative, 

„, Lxd HL 2 -L(L H) 1 fdL \ . 

M = — + aD - + ^U x£ )' (30) 

where L — \L\. Within the sigma-model approximation, the length of the vector L should 
be treated as a constant, L 2 X + L 2 + L 2 = const ~ (2M ) 2 . Thus, in the linear approximation 
L x ~ 2M — (ly+ll) /4M and the two components, l y and l z can be considered as independent 
variables. It is in line with our experimental observation that the component l x is completely 
negligible. The effective anisotropy energy can be taken in the form 

1 < w2 1 f d2 
2(93-^+2(^4 

where the additive constant is omitted. Free oscillations of the two components at H = 
correspond to two independent magnon modes (F- and AF- modes), described by the 
following equations 

d 2 l z 2 A d dl z 

+ u F i z = 0, m = x— L z + y- 



W.(X) = - (q 3 - qi ) l 2 z + ^[--Qi)il (31) 



dt 2 * ' A 2 7 AM dt 

d 2 l y 2 , 1 dl 



dj2+w | pis = 0,m = -z— J. (32) 

Now let us discuss the excitations of the modes by light pulses. The interaction of 
the spin system with the light is described by the Hamiltonian ffTO]) . that for the specific 
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case of circularly or linearly polarized light reads as ( TTTj) or (Tl2|) . respectively. Within 
the sigma-model approach, for different polarizations the interaction terms enter different 
parts of the Lagrangian (1291) : the circularly polarized light contributes to the effective field 
H = H e ' a , whereas the effect of the linearly polarized light is described by the time- 
dependent contribution 

SW a (L,t) = ^(*)£*(*)&y W L fc L,, (33) 

to the effective anisotropy energy W a {L). Among all these contributions to the Lagrangian, 
we need to find terms linear on l y and l z , which produce the "driving force", i.e., lead to a 
non-zero right-hand side in the equations of motion (132]) . 

The light-induced effective field is directed along z-axis, and it is easy to see that the 
term (H ■ (d x L)) gives no "driving force" contributions for both modes. The gyroscopic 
term with dl/dt provides such a term for y-component of the vector Z, proportional to 
if| ff ' CT± L x (dl y /dt), but not for its z-component. Thus, for the state of interest (L = L x x in 
the ground state), the IFE can excite the AF-mode only. In the discussion presented above, 
the only part proportional to b xyxy L x l y gives an essential contribution to 5W a (L, t). Using 
these relations, one can find that all the terms do not affect the equation for l z (F-mode), 
whereas the equation for l y describing the AF-mode acquires nonzero right-hand side and 
reads as 

+ "l F ly = -1lM,-^- + A(2 1 M )X 1 , (34) 

where h y n = —d5W a (L,t)/dl y is the effective field. Then, after the delta function substitu- 
tion S(t)S*(t) — > IoS(t), we arrive at the following initial conditions for this equation 



lM f xyz I° 

\ L y)t=o ~ ■ J - 

dl 



(ly) t=0 = ± ' °T ° , (35) 



2 

-4 7 2 AM*& w Csin20, (36) 



dt 



t=o 

where Jq and determine independent action of circularly and linearly polarized light, 
respectively, with Iq ± and Iq 11 being the corresponding integrated pulse intensities. As one 
can see from the equation, within the sigma-model approach the effective magnetic field 
created by the IFE enters the equation through its time derivative only, whereas the inertial 
mechanism is caused solely by ICME. Thus the field-derivative mechanism of the action 
of IFE, discussed previously for compensated antiferromagnets,- 1 ^ is responsible for the 
excitation of spin oscillations in the r 4 -phase of dysprosium orthoferrite investigated here. 
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We conclude that it is difficult to realize the inertial mechanism of the field pulse action 
in the majority of orthoferrites at high temperatures where the same I^-phase is present. 
The inertial mechanism has been observed for a special phase of holmium orthoferrite where 
the vector L is not collinear with the symmetry axis.— On the other hand, for the present 
experiment the ICME leads to inertial mechanism of the spin excitations. 

After the action of the pulse, only free spin oscillations persist in the system. They are 
described by the solution 

l% ee = acos(u AF t + £), (37) 



where the amplitude a and the phase £ are determined by the initial conditions (1361) as 
follows: 

4AM b xvxv I^ sin 29 . . 

tan£ = ^ ° l yx y . 39 

f xy JS \/d 2 -Aq x 

a -w,./^-- 7 - ) 2 | (2^M ^sin2^ 
V 4 d 2 — Aqi 

Finally, we arrive at the previous result: if one of the two mechanisms, IFE or ICME, 
is dominating, the the phase of the m z oscillations takes the values £ = 0, it or £ = ±7r/2, 
respectively. Thus, the observed time dependence of the Faraday rotation oscillations is 
proportional to sinc^AF^ or coswaf^ for the dominating role of IFE or ICME, respectively. 
If none of the mechanisms is truly dominating, then the observed phase should take an 
intermediate value given by Eq. ( 1401) . 

It is worth to note that the condition for domination of a certain effect does not translate 
into a plain comparison of the effective constant values f xyz and 2M b xyxy for IFE and ICME, 
respectively. The point is, the ICME contributes through the inertial mechanism that is 
much more effective than the field-derivative mechanism involved in the action of IFE. In 



our calculation, this leads to appearance of the large multiplier A/ y ' d 2 — Aq\ = jH ex /uAF, 
where 7-ff ex ~ 20 THz, H ex = 2AM ~ 600 T is the exchange field of orthoferrite,— in the 
contribution of ICME, see Eq. ( 1401 . Therefore the domination of IFE, for the same value of 
the pulse fluence, needs at least 50 times higher value of the corresponding constant, and 
the ratio f X yz/2M b xyX y is expected to be large enough for orthoferrites. Thus, the above 
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analysis gives us a possibility to estimate the values of constants responsible for different 
inverse magneto-optical effects, IFE and ICME. 

C. Comparison between the theory and the experiment 

In the previous discussion, based on the Landau-Lifshitz equations and the nonlinear 
sigma model, we came to the conclusion that the time dependence of m z induced via IFE 
and ICME is proportional to sin waf^ and cos waf^, respectively. The phase of the oscillation 
is constant and m z is proportional to either sin c^af^ o r cos waf^ i n some region of the pump 
wavelength in Fig. 0(a). When the pump pulse is in the visible region (<800 nm), the probe 
polarization and m z oscillate as sinajAF^- This property is independent of temperature as 
shown in Fig. [8] (b). On the other hand, when the pump pulse is in the near- infrared region 
(1000-1100 nm), the probe polarization and m z oscillate as cosojaf^- Thus, we can conclude 
that the visible and near-infrared light pulses dominantly induce spin precession via IFE 
and ICME, respectively. 

A number of reasons can be given for why the dominant effect varies with pump wave- 
length. IFE is induced by a pulse whose wavelength is near the transition 6 A 1 ->■ 4 T 2 at 700 
nm. On the other hand, ICME is induced by a pulse whose wavelength is near the transition 
6 Ai — > 4 Ti at 1000 nm. In addition, the Faraday rotation angle increases with decreasing 
wavelength in DyFeOs.— ^ This tendency agrees with the result for the IFE. 

V. CONCLUSIONS 

We have studied the dependence of photoinduced spin precession in DyFeOs °n the 
wavelength and polarization of a pump pulse with a pump-probe magneto-optical technique. 
The polarization rotation of the probe pulse was dependent on the pump polarization. Pulses 
propagating along the z-axis with both circular and linear polarizations induced an effective 
magnetic field (IFE and ICME) and spin precession. The dominant component of the 
dielectric permittivity in both effects was e xy , and IFE and ICME were induced by its 
antisymmetric and symmetric parts e a xy and e s xy , respectively. 

The phase and amplitude of the spin precession were dependent on the pump wavelength 
in DyFeOs. A difference in the pump wavelength changes the dominant effect, giving rise 
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to the spin precession. A visible pulse (wavelength <800 nm) induced the IFE, and the 
oscillation of the probe polarization was proportional to sinwAF^- On the other hand, a 
near- infrared pulse (wavelength of 1000-1100 nm) induced the ICME dominantly, and the 
oscillation was proportional to cos waf^- When the pump wavelength was near the electron 
transition & A l -> 4 T 2 at 700 nm and 6 Ai — > 4 T X at 1000 nm, the amplitude of the oscillation 
was larger than that of the other region. 

The ratio of the effective magnetic fields via IFE and ICME, f X y Z /2Mob xyX y, is expected 
to be large enough for orthoferrites. However, the ellipticity of spin precession with AF- 
mode is also so large. Therefore, even though linearly polarized light pulse induces so weaker 
magnetic field than circularly polarized one, ICME can give the same order contribution as 
IFE. 
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FIG. 1. (a) The static magnetic structure of DyFe03, with the four Fe 3+ spins regarded as sat- 
isfying Mi ~ M% and . -antiferromagnetic spin 
precession. 
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FIG. 2. The experimental setup geometry. BBO was used for frequency doubling of the pump 
pulse, if necessary. WP: Wollaston prism, GTP: Glan-Taylor prism, P: polarizer, A/2: half-wave 
plate, A/4: quarter- wave plate. 
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FIG. 3. Pulse polarizations. Circularly polarized pulses are denoted . Linearly polarized pulses, 
denoted LI, L2, L3, L4, L5, and L6, were tilted at — vr/4, 7r/4, 0, 7r/2, and =pa from the x-axis, 
respectively, where tana = 2. 
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FIG. 5. Oscillation of the probe polarization 6{t) as a function of the time delay between the pump 
and probe pulses. Six types of probe polarizations were used: (a) circular polarization a^, and 
linear polarizations tilted at (b) =F45°, (c) 0°, and (d) 90° with respect to the x-axis. 
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FIG. 6. Time-resolved pump-induced probe polarization 0(t). The pump polarization was circular, 
and the graph shows the shift of the probe polarization when changing the pump polarization from 
cr + to a~ . The pump wavelengths were (a) 750 nm and (b) 1050 nm, and the probe polarizations 
were L5 and L6. 
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FIG. 7. (a) Initial phase and (b) amplitude of the oscillation of the polarization as a function of 
the pump wavelength. The amplitude is normalized by the pump fluence. 
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FIG. 8. (a) Temperature dependence of the induced spin precession frequency. Pump wavelengths 



were 750 nm and 1200 nm. The solid line shows the magnon frequency taken from Ref. 
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(b) 



Temperature dependence of the initial spin precession phase for the pump wavelength of 750 nm. 
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